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A colloidal system of spheres interacting with both a deep and narrow attractive potential and a shallow 
long-ranged barrier exhibits a prepeak in the static structure factor. This peak can be related to an additional 
mesoscopic length scale of clusters and/or voids in the system. Simulation studies of this system have revealed 
that it vitrifies upon increasing the attraction into a gel-like solid at intermediate densities. The dynamics at 
the mesoscopic length scale corresponding to the prepeak represents the slowest mode in the system. Using 
mode coupling theory with all input directly taken from simulations, we reveal the mechanism for glassy arrest 
in the system at 40% packing fraction. The effects of the low-g peak and of polydispersity are considered in 
detail. We demonstrate that the local formation of physical bonds is the process whose slowing down causes 
arrest. It remains largely unaffected by the large-scale heterogeneities, and sets the clock for the slow cluster 
mode. Results from mode-coupling theory without adjustable parameters agree semi-quantitatively with the 
local density correlators but overestimate the lifetime of the mesoscopic structure (voids). 

PACS numbers: 82.70.Dd 61.20.Ja 64.70.Pf 



I. INTRODUCTION 

Many materials solidify into nonequilibrium structures 
when particle interactions become strong and if crystalliza- 
tion proceeds too slowly. Colloidal dispersions provide a 
number of model systems to study these still little under- 
stood solidification phenomena as particle potentials can be 
tailored and detailed experimental observations are possible. 
For example colloidal hard spheres have been studied inten- 
sively, where addition of nonadsorbing polymer induces an 
attraction, whose strength and range can be controlled d. 
The case of short-ranged attractions has turned out especially 
rich. Colloidal particles interacting with already moderately 
strong short-ranged attractions can form metastable amor- 
phous solids which exhibit exceedingly long life times. 

Depending on the colloidal density and attraction strength 
different types of metastable arrested solids can be formed. 
At high densities glasses are observed when repulsions hinder 
and even prevent structural rearrangements |2fl. Weak short 
ranged attractions at first melt these 'repulsion-driven glasses' 
because they distort and loosen the local packing. At attrac- 
tion strengths somewhat higher, yet still of the order of the 
thermal energy, physical bonds are formed in dense systems, 
which leads to aggregation into 'attraction driven glasses' 
JH IU, Ht]. At attraction strengths high compared to thermal 
fluctuations aggregation phenomena proceed far from equi- 
librium at low density, resulting in tenuous solids, i.e. gels. 
At even lower density, the aggregation process leads to cluster 
formation and aggregation 00111O- 

The connection between repulsion and attraction driven 
glass transitions at high densities has been understood within 
a microscopic theoretical framework, namely mode coupling 
theory (MCT) HI [H [H [H GJ]. Yet, the mechanisms 
of solidification at intermediate attraction strengths and low 
to intermediate densities are still not completely understood 



lfl6i rTvh - Especially the interplay of equilibrium phase tran- 
sitions and a ggre gation effects is presently under scrutiny 
0, HI EH IT9L l20ll . Another aspect that has been considered 
at low densities is the existence of weak long-ranged repul- 
sive interactions ('barriers') between the particles, that may 
induce density modulations (or non-percolating clusters at low 
density), and thus prevent solidification for some parameter 
ranges lEltl . Small wavevector scattering peaks indicative of 
mesoscale correlations have also been observed and employed 
to characterize the cluster or gel structures [fToL I22I l23ll . Their 
observation in equilibrium systems with long-ranged barriers 
is expected because these systems show tendencies to micro- 
phase separate and to form interconnected structures of clus- 
ters and voids J3- Yet, such long distance correlations and 
clusters of finite size were also observed in systems where 
barriers supposedly are absent and there presumably have a 
nonequilibrium origin and arise from phase separation and/or 
aggregation and coarsening JHlit]. Nonequilibrium origins of 
the low angle scattering peaks also were suggested by the ob- 
servation of their time evolution 1I22I l25ll . Alternatively, to 
reach gelation from the fluid, liquid-gas phase separation can 
be prevented setting a maximum number of bonds per par- 
ticle ll26ll or total number of bonds in the system ll27ll . In 
these cases, gelation is directly connected to percolation, and 
the low-q modes facilitate the relaxation of the whole system, 
due to the lack of stiffness (only two to four neighbours are 
allowed on average). At higher density, such modes are ex- 
pected to be less important than collective relaxations, and 
a general theory should account for both relaxation mecha- 
nisms. 

In this contribution we test quantitatively the MCT predic- 
tions for low density attractive glasses or high density gels 
and study the role of the cluster or void structure on bond for- 
mation at an intermediate density. The input structure factor 
needed for the MCT calculations is taken from simulations, 
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and we compare the properties of the density correlation func- 
tions. Additionally, by comparing the results of MCT calcula- 
tions for systems with and without small wavevector pre-peak 
in the structure, we highlight the importance of mesoscopic 
heterogeneities on attraction driven dynamic arrest. We con- 
sider a system of particles interacting with a narrow attraction 
and a weak long-ranged repulsion whose dynamics has been 
studied intensely by simulations yl|2ja|2ja|30[]. At the consid- 
ered density of 40% packing fraction, our system is above the 
percolation threshold, and exhibits an equilibrium structural 
pre-peak at small wavevectors for parameter ranges where the 
barrier suppresses phase-separation. By switching off the re- 
pulsive barrier, the cluster or void structure can be eliminated 
and the system becomes homogeneous (below phase separa- 
tion). 

Similar studies, taking the S q from simulations as input for 
MCT, have been performed previously in repulsion driven sys- 
tems, for instance Lennard-Jones [f3lL l32ll . molecular glasses 
ll33tl . silica ll34tl or sodium silicate melts l35ll . For attraction 
driven colloidal glasses, a similar study has recently been per- 
formed by Manley et al. using experimental structure fac- 
tors lfl8ll . and we reach an analogous conclusion concerning 
the mechanism of dynamic arrest. Our study goes beyond 
Ref. Il8ll because the input to our MCT calculations is taken 
from simulations of the system we address without adjustable 
parameters, while Manley et al. adjusted the colloid density 
arbitrarily and assumed a specific form for the static struc- 
ture factor at large wavevectors. We find that it is exactly the 
structure at large wavevectors that dominates bond formation. 
Similar effects in the dynamics related to a prepeak in S q have 
also been observed in sodium silicate melts ll35i[36ll . 

II. MODE COUPLING THEORY 

This section aims to give a brief overview of the most 
important results within MCT concerning the description of 
liquid-to-glass nonergodicity transitions l37l l38il . We focus 
on the coherent part of the density correlation function as it 
provides most insights into the physical mechanisms causing 
glassy arrest. 

MCT gives a self-consistent equation of motion of the (nor- 
malized) intermediate scattering function &q(t), which is the 
coherent part of the autocorrelation function of density fluctu- 
ations with wavevector qof N particles, defined by 

*«( f ) = IJT t <exp{/q[r ; (r)-r 7 (0)]}). (1) 

JVZ> q j j=1 

The normalization to unity at time t = is provided by the 
static structure factor S q = £f ; = i(exp{iq[r,- — rj]})/N, which 
captures equilibrium density correlations. The equation of 
motion for ^ q (t) takes the form of a relaxation equation, 
where retardation effects with respect to exponential relax- 
ation on diffusive timescale T q = ( y D s q 2 /S q ) 1 are contained 
in a memory kernel m q (t), 

<T,* 9 (f )+*,(*)+ f dt'm q {t-t ')% (f')=0 (2) 

J 



The initial decay constant 1 q describes short time diffusive 
particle motion, and is set by the short time collective diffu- 
sion coefficient D s ; it captures instantaneous particle interac- 
tions and will not play an important role at the glass transition. 

The central quantity capturing slow structural rearrange- 
ments close to glassy arrest is the memory function m q (t) 
which is given in MCT-approximation: 

m q (t) = *•,[{*(*)}] = i^Vqkp **(*)*p(0 0) 

L kp 

The vertices Vqk p couple density fluctuations of different 
wavelengths and thereby capture a nonlinear feedback mecha- 
nism in dense fluids, which is interpreted as 'cage effect' l38ll . 
The memory kernel can be regarded as a generalized friction 
kernel, as can easily be verified after time-Fourier transfor- 
mation of Eq. (0. MCT is a first principles approach as the 
vertices are calculated from the microscopic interactions 

p 2 

V qk p = S q S k S p -^ [q • kc* + q pc p ] 2 8(q k p) (4) 

The mode coupling approximation for m q yields a set of equa- 
tions that is solved self-consistently. Hereby the only input to 
the theory is the static equilibrium structure factor Sq, which 
enters the memory kernel m q directly and via the direct corre- 
lation function c q = (1 — l/S q )/p, with p = N/V the average 
density. According to MCT the dynamics of the dense fluid 
is therefore completely determined by equilibrium quantities 
plus one scale factor to set the time scale. In the following, we 
will take S q directly from simulations. 

Solutions of Eqs. (04]) show a bifurcation scenario due to 
the nonlinear nature of the equations. The bifurcation point is 
identified with an idealized liquid-to-glass transition. A quan- 
tity of special interest is the long-time limit of the normal- 
ized density correlator, f q = lim f _»,x, < I > , ? (f), often referred to 
as glass form factor or Edwards-Anderson nonergodicity pa- 
rameter. It describes the frozen-in structure of the glass and 
obeys 

TZT = Ti\S\- (5) 
1 Jq 

In the fluid regime density fluctuations at different times 
decorrelate, so that the long time limit vanishes, f q = 0. On 
approaching a critical packing fraction § c or a critical tem- 
perature T c , MCT states that strongly coordinated movements 
are necessary for structural rearrangements to relax to equilib- 
rium. MCT identifies two slow structural processes, p- and a- 
process, when the glassy structure (described by f q ) becomes 
metastable and takes a long time to relax. In the idealized pic- 
ture of MCT, dynamical arrest sets in at the glass transition 
with f q > when the particles are hindered to escape from 
their neighbouring environments. This also is accompanied 
by diverging relaxation times. Glassy states are called noner- 
godic states in MCT. The value of the glass form factor at the 
transition is called critical nonergodicity parameter . 

Although experiments on molecular glass formers have re- 
vealed that the dynamics very close to the transition point is 
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dominated by thermally activated hopping processes, which 
the described (idealized) MCT cannot account for, MCT has 
been very successful in describing the approach to glassy ar- 
rest. It gives a quite accurate description of structural re- 
laxation in colloidal dynamics and there especially of the a- 
process. 

For liquid states and large times the correlator approaches 
the a-scaling law <£> g (t) — > 3> 9 (f) and f = f/x. Here the scal- 
ing functions <P q (t) are independent of temperature or other 
control parameters. The oc-relaxation scale x is given in MCT 
by 

x = x(e) =x /|e| Y withy = -^- + ^- (6) 

2a 2b 

and depends only on the separation e = Ssh- from the criti- 
cal point. The scaling factor Xo needs to be determined from 
matching the microscopic dynamics and the p-scaling regime 
to the a-scaling regime. The anomalous exponents a and b are 
solutions of the equation 

r 2 (l-jt-)/T(l-2x) = X (7) 

for < x = a < 1/2 and — 1 < x = —b < 0. The exponent 
parameter X enters in the second order of the asymptotic ex- 
pansion of the right-hand side of Eqs. ( 1315b around the critical 
plateau f q . Therefore it depends on the structure factor S q at 
the critical point via the vertex Eq. (@}. 

In the vicinity of the critical point, von Schweidler's power- 
law describes the initial a-relaxation from the plateau to zero. 
It is nothing more than the short-time expansion of the oc- 
master functions, which is up to second order 

%(t) = f q - h q t\\ + k q t b ) + o (P). (8) 

The coefficients h q are called critical amplitude, the k q are 
correction amplitudes rf39tl . Von Schweidler's law is the origin 
of stretching (viz. non-exponentiality) in the a-process of 
MCT. 

The final decay of the structural relaxations in different cor- 
relators provides a definition of the a-relaxation times. We use 
®q(t = T ? ) = fq/20. In MCT the increase of the relaxation 
times in different correlators is strongly coupled. Their di- 
vergence is inherited directly from the diverging a-relaxation 
scale x of Eq. (|S). For the a-relaxation times of the density 
fluctuations in & q (t) this means a separation into a-scale and 
a dimensionless factor t q containing the wavevector depen- 
dence. 

z q =i q x(e). (9) 

For large wavevectors, L ~ q~ l l h holds in MCT ll40ll . 

Summarizing this short presentation of MCT, let us note 
that the wavevector dependences of the various amplitudes in 
the asymptotic MCT predictions will enable us in the follow- 
ing to identify the physical mechanisms causing glassy arrest. 
More details about MCT, the asymptotic expansions and the 
scaling-laws can be found in 1137t I38L I39L 14111 . 



HI. SIMULATION SETUP 

Molecular dynamics simulations in the canonical ensemble 
were performed considering 1000 quasi-hard particles inter- 
acting by a short range attraction. Because we aim to study the 
fluid to non-ergodic transition induced by attractions, equilib- 
rium phase transitions, i.e. crystallization and liquid-gas sepa- 
ration, were suppressed by suitable choices for the interaction 
potential, which we introduce in the following. 

The short range attraction mimics the interaction between 
colloidal particles induced by non-adsorbing polymers in a 
colloid polymer mixture. For monodisperse colloids, this 
attraction is given by the Asakura-Oosawa interaction Ql, 
which is slightly modified to include polydispersity ll42ll . The 
attraction strength is set by the concentration of polymers, § p , 
and the range by the polymer size, t, (see below). This po- 
tential has been slightly corrected near contact r = dn to en- 
sure that the total interaction potential has the minimum at dyi 
(d\2 = (a\ +02), with a\ and ai the radii of the interacting 
particles) I29I 
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FIG. 1: Interaction potential for two particles with the average ra- 
dius. The polymer fractions are ty p = 0.42 and ty p = 0.25 for the thin 
and thick curves, respectively (values close to the glass transition in 
the simulation and MCT, respectively). Note that in our units, the 
thermal energy is IcrT = 4/3. 

Crystallization is avoided by using a polydisperse system: 
particle sizes are distributed according to a flat distribution of 
width 8 = 0.1a, where a is the mean diameter. The core-core 
repulsion is given by V sc (r) = kgT (r/dn) 36 - 

At high polymer fractions, or attraction strength, this sys- 
tem separates in two fluid phases with different densities, di- 
lute and dense - the critical point is at (j)p P w 0.29. To avoid 
this transition, which would interfere with the attractive glass, 
a repulsive long-range barrier has been added to the total in- 
teraction potential, which extends from r = c/12 + 2c; to r = 4a, 
and its height is only IkgT (equal to the attraction strength at 
§ p = 0.0625). The barrier raises the energy of a dense phase, 
so that liquid gas separation does not take place. The resulting 
total interaction potential, V to t = V sc + Vao + Vbar> is analyti- 
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cal everywhere and allows straightforward integration of the 
equations of motions. The total interaction potential is pre- 
sented in Fig. [T]for particles with the average radius. 

The inclusion of the repulsive barrier in the interaction 
potential effectively inhibits liquid-gas separation ll29tl . but 
causes holes and tunnels in the system. This is reflected in 
the structure factor as a low-q peak, which grows and moves 
to lower-g values, as the strength of the attraction (namely, § p ) 
increases. The effect of this barrier on the glass transition is 
studied below within MCT 

Lengths are measured in units of the average radius, a, 
time in units of y / 4a 2 /3v 2 , where the thermal velocity v was 
set to \/4/3. Equations of motion were integrated using the 
velocity-Verlet algorithm, in the canonical ensemble (constant 
NTV), to mimic the colloidal dynamics. Every n, time steps, 
the velocity of the particles was re-scaled to assure constant 
temperature. No effect of n t was observed for well equili- 
brated samples. The time step was set to 0.0025. The range of 
the attraction, 2t,, is set to 2% = 0.2a. The density of colloids 

is reported as volume fraction, <j) c = ^fta 3 ^1 + (j^j ^ n c , 

with n c the colloid number density, and the attraction strength 
is given by the polymer volume fraction, p (with i; = 0.1, 
the minimum of the attraction for average sized particles, at 
r = 2a, is V mtn = -16k B Tty p ). 

The dynamics of this system has been analyzed previously 
within the framework of MCT, i.e. using the density corre- 
lation functions |29[ l43ll . Increasing the attraction strength, 
§ p , at <j) c = 0.40 a glass transition is obtained at § c = 0.4265, 
which shows the qualitative features of attractive glasses, as 
predicted by MCT. The critical parameters given below for 
the transition were obtained analysing <& q (t) for the fluid 
state (j) p = 0.42, close enough to the glass to show the typi- 
cal two-step decay (see Fig. [H) - 1500 independent config- 
urations were used to calculate &q(t), from 15 independent 
"quenches" from hard spheres at the same density. The p- 
regime and early a-decay were analyzed from t = 2 to t = 500. 
A wavevector-independent von Schweidler exponent was fit- 
ted using the correlators at all wavevectors, whereas the non- 
ergodicity parameter, f c q and critical amplitudes h q and k q 
were actually fitted for every wavevector. 

The structure factors needed as inputs to MCT were cal- 
culated from the definition of S q , using only the allowed 
wavevectors q = 2Tt/L(n x ,n y ,n z ), with L the box size and n x , 
n y and n z integers. Starting from q = 2u/L, the next value 
of the q modulus is selected if the ^-separation is larger than 
O.lfl -1 , up to qa = 40. The structure factors were then inter- 
polated to have a constant g-grid. 

IV. ASPECTS OF THE NUMERICAL MCT SOLUTIONS 

For the numerical solution of the MCT equations, algo- 
rithms were used that have been developed in the recent years 
& 

Dynamic and static analyzes were performed by itera- 
tively solving Eq. (O with the memory functional given by 
Eqs. ( 1314b . The results were accepted if a convergence to 



<£>q(t) with a relative accuracy of 10~ 15 was achieved at each 
t,q. To extend the calculation onto logarithmic time scales 
without running into inefficient time-discretization for late 
times, we used an algorithm for the convolution in Eq. (0 
that doubled the initial time step of 10~ 9 D 2 /D s every 256 time 
steps, where D = 2a is the particle diameter. The critical poly- 
mer concentrations were attained by a bisection method and 
determined up to a relative accuracy of 10~ 5 in concentration. 

The structure factors for the MCT calculation input were 
taken directly from the simulation and linearly interpolated 
without further smoothing. We used a wave vector grid with 
M = 400 grid points and a cutoff of qD = 80. From it is 
known that with these values neither the discretization nor the 
cutoff influence the results significantly. For very small wave 
vectors qD < 0.3 the algorithm does not produce the correct 
results because of numerical error propagation. This can be 
recognized by the static results shown in Figs. [4]|5] for the 
nonergodicity parameter f q , in Fig. [6] for the critical ampli- 
tudes h q and in Fig. |7]for the a-relaxation times. The plotted 
results exhibit a sudden drop down to for q — ► 0, whereas 
they should take a finite value. Despite this, the results for 
larger q are not invalidated, as was verified by removing the 
incorrect values from the integrations in Eq. (|3). 



V. RESULTS AND DISCUSSION 

In order to quantitatively describe the results from the sim- 
ulations for the polydisperse system of particles interacting 
with the potential shown in Fig. Q] three different calculations 
within MCT were performed. They differ in the input static 
structure factors, all of which were obtained directly from 
simulations. Three additional simulation studies were per- 
formed solely to generate the S q where the pair potentials em- 
ployed in the different simulations differed. Thus, we could 
highlight the importance of (/) particle polydispersity, and (if) 
the long-ranged repulsive barrier, in the structural relaxation. 

For convenience and clarity we name the different MCT 
calculations in the following way: System (A) is the monodis- 
perse model without repulsive barrier, whereas we refer with 

(B) and (C) to the monodisperse and polydisperse systems, 
respectively, with repulsive barrier. While the calculation in 

(C) thus uses exactly the S q of the system whose dynamics 
we aim to describe, the MCT we use considers a monodis- 
perse system. Calculation (C) thus also is only approximative. 
True multi-species MCT calculations like in Ref. [46] would 
be required to capture all polydispersity aspects, yet are too 
demanding in the present case. All elements of the matrix 
of partial structure factors would be required, and the multi- 
ple wavevector integrations over the required large ^-ranges 
would crucially slow down the MCT numerics. In order to 
stress the approximative character of the MCT calculation (C), 
we call the simulation, where the dynamics is analyzed, sys- 
tem (D). Some results of the calculations like critical polymer 
concentrations, exponent parameters, exponents, and localiza- 
tion length (in units of the diameter D) are summarized in ta- 
ble u 
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(A) MCT 





No 


0.2356 


0.752 


0.555 


0.00501 


(B)MCT 





Yes 


0.2464 


0.759 


0.544 


0.00510 


(C) MCT 


0.10 


Yes 


0.3646 


0.775 


0.517 


0.02628 


(D) Simulation 


0.10 


Yes 


0.4265 


0.863 


0.37 


0.0158 



TABLE I: Critical polymer concentrations, exponent parameters, 
von-Schweidler exponents and localization lengths resulting from 
the different MCT computations: (A) monodisperse without barrier, 
(B) monodisperse with barrier, (C) polydisperse with barrier. The 
first columns give polydispersity 8 and whether a repulsive barrier 
exists. Listed under (D) are the corresponding parameters from the 
analysis to the polydisperse simulation with barrier of the fluid state 
with § p = 0.42. All states are at colloid packing fraction <j) c = 0.40. 



A. Structure factors 

To clarify the influences of the different pair potentials and 
the polydispersity on the equilibrium structure we compare on 
the one hand structure factors for the same parameters ((j) c -, § p ) 
for the different potentials and on the other hand results at the 
MCT-critical points that are obtained for the different poten- 
tials. The Sq in Fig. [2] for § c = 0.40 and = 0.25 show all 
a primary peak at qD ml. 5 that indicates the local fluid or- 
der. A peculiarity, which is often seen in gels of intermediate 
and lower densities, is a low-q peak in the static structure fac- 
tor. It appears on the length scale of the voids in the structure 
when the sol-gel transition line is reached, increases towards 
the transition flS^l, and shifts to slightly lower q- values (cf. 
Fig. [3}. A common interpretation of the low-g peak is that 
it indicates the onset of an arrested phase separation at higher 
attraction strengths. At lower density, below the percolation 
threshold, this peak marks the presence of clusters in the sys- 
tem, although similar internal structures for gels and indepen- 
dent clusters have been reported ifioll . The systems (B) and 
(C) with repulsive barrier show this low-g peak, whereas it is 
absent in the model (A) without repulsive barrier. The latter S q 
clearly grows in the limit q — > for increasing polymer frac- 
tions, indicating the proximity of the liquid-gas critical point. 
The S q with barrier stay finite as a result of the weak repulsion 
preventing phase separation; rather they develop the prepeak 
which may signal closeness to microphase separation. 

Besides the long-ranged barrier, polydispersity has impor- 
tant effects on the equilibrium structure. The inset of Fig. [2] 
displays S q from qD = 30 — 80 for the different systems. The 
Sq of the monodisperse systems both with and without bar- 
rier virtually coincide and show distinct oscillations for these 
wave vectors, unlike in the polydisperse case. Indeed Sq for 
(C) starts to deviate from (B) above qD w 12 and decays to the 
noise level above qD w 45. This rapid decay to unity in the 
Sq of the polydisperse systems is due to slight differences in 
the distances where the (partial) pair correlation functions for 
differently sized particles show their contact maxima. This 
distribution in the contact distances leads to negative inter- 
ferences in the oscillatory large-g pattern in Sq, which gets 
canceled in the averaged structure factor of the polydisperse 
system. This effect of the short-ranged attraction, viz. the in- 
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FIG. 2: Sq from MD-simulations at a colloid packing fraction of 
§ c = 0.40 and a polymer concentration of § p = 0.25, which is in the 
gel close to the critical point in MCT. The repulsive barrier affects 
Sq in the region below qD < 13; system (A) without barrier (dashed- 
dotted, black) shows neither a prepeak nor a primary peak which is 
as high as in systems (B,C) with barrier (full red, dashed black). The 
inset demonstrates that polydispersity causes the q-tail oscillations 
in S q to be suppressed for large wavevectors. The S q for the polydis- 
perse model (C) falls below the noise level for qD > 45, whereas the 
monodisperse systems (A) and (B) virtually coincide there. 



creased probability of particle contact, is thus only contained 
in the S q of systems (A) and (B). 



♦ ■ -» (D) simulation 




FIG. 3: Critical Sq at the boundaries of the gel phase: Black dashed- 
dotted, red full and black dashed lines mark results from (A), (B) 
and (C). The red diamonds indicate S q at fy p = 0.42, which is close 
to the arrested state in the simulation (D). In the simulation, the pre- 
peak at qD ~ 2.5 rises and shifts to smaller (/-values with increasing 
polymer concentration. The inset gives an enlarged view of the q- 
tail where systems (A,B) show pronounced oscillations driven by the 
short-ranged attraction. 

Fig. [3] shows Sq at the critical points in MCT (see table |I| 
and in the simulation ((D), ty p ps 0.4265). It was possible to 
take the structure factors for the MCT input directly from the 
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monodisperse simulations, because MCT predicts the states 
to be nonergodic already at rather low §„, while the actual 
system only freezes at higher attraction strengths. The fea- 
tures discussed in respect to Fig. [2] can be recognized again 
even though the attraction strenghts vary. While S q of the sys- 
tem (A) without barrier differs from the other systems at small 
wavevectors, polydispersity forces the averaged S q of systems 
(C,D) to approach unity quickly at large wavevectors. 



B. Nonergodicity parameters 

From the equilibrium structure factor input of the differ- 
ent systems we calculated the critical Edwards-Anderson non- 
ergodicity parameters f q using Eq. (f5]). The bifurcation oc- 
curs for the colloid packing fractions of (j) c = 0.40, which has 
been used throughout the analysis, at polymer concentrations 
if p = 0.2356 (A), ty c p = 0.2464 (B) and <j>£ = 0.3646 (C) re- 
spectively. Note that the critical polymer concentrations are 
lower than the critical tfL in the simulation, system (D). Luck- 
ily, the transition point in system (A) occurs below the criti- 
cal point, which suppresses effects from the liquid-gas sepa- 
ration. For system (B), crystallization is far too slow to affect 
the results (although this state is indeed metastable with re- 
spect to crystallization). The trend of MCT to overestimate 
the tendency to freeze leads almost to a factor of 2 in terms of 
attraction strengths l45ll . This discrepancy has also been ob- 
served in former comparisons of MCT with binary Lennard- 
Jones fluids l3~lll . 

The non-ergodicity parameter, f q , basically oscillates in 
phase with the S q . The shape of f q serves to identify the 
leading mechanism for the freezing. A repulsion driven tran- 
sition creates an f q with pronounced peaks and lower val- 
ues for small wave vectors. Characteristically, it decreases 
quickly to zero for increasing q. The width of f q as function 
of q can be taken as a measure for the localization length that 
describes the spatial extent that a single particle can explore 
within its glass cage. For repulsion driven glass transitions 
one generally finds a localization length of the order of the 
Lindemann-length, viz. a value around a tenth of the aver- 
age particle separation l38ll . If attraction drives the transition, 
the critical nonergodicity parameters have higher values and 
smaller localization lengths, showing up in much wider f q - 
distributions than for repulsion driven nonergodic states. The 
width of f q as function of q is now set by the attraction range. 
In our analysis f q compares quite well with the simulation re- 
sults for those models where an attraction driven large-q tail 
is present in S q (see Fig. |4|. The MCT calculation for the 
polydispersity-smeared-out S q gives far too small f q and too 
large localization lengths what resembles more results from 
repulsion driven systems. We conclude from the agreement 
between the fjj from simulation and the MCT calculations 
with attraction-driven large-g tails in S q , that the simulations 
exhibit an attraction driven glass transition or what could be 
referred to as a gel. Note that this good agreement is found 
despite the overestimate of the critical attraction, i.e. the tran- 
sition is wrongly located but its principal property is correctly 
predicted. 




' 20 40 60 80 

qD 

FIG. 4: Critical nonergodicity parameters f q at the transition in the 
simulation (D) (red diamonds), from MCT with monodisperse (A,B) 
S q (black dashed-dotted and red full line) and polydisperse (C) S q 
(black dashed line). Serious differences occur, when the average S q 
of a polydisperse system is used as input to a monodisperse theory; 
the shape of f q for (C) resembles more the one of a repulsive glass. 



An important point to be checked is the role of the long 
ranged repulsive barrier. We have already stated above that 
there is hardly any difference between the critical polymer 
concentrations in the monodisperse calculations without (A) 
and with barrier (B). Fig. [5] highlights the influence of the 
barrier on the nonergodicity parameters f q . The localization 
length as well as the attraction driven character of the glass re- 
main unchanged. A significant effect could only be observed 
in the vicinity of the prepeak and the primary peak, where S q 
also changes. Domains of higher wave vectors stay practically 
unaffected (see also Fig.|4|. This indicates that the modes on 
the low angle q-peak, which is related to the void structure 
seen in the simulations, follow the relaxation of an attraction 
driven glass without dominating it. It is dominated solely by 
large-g modes. 



C. Critical Amplitudes and Von-Schweidler-Law 

Asymptotic expansions of Eq. (O around the critical 
plateau f c q introduce the critical amplitude h q in linear order 
of the so-called p-process. During the p-process, the dynam- 
ics on all length scales is strongly coupled. The amplitude h q 
measures the participation of the correlator at wavevector q in 
this process. The p-process describes the rate limiting process 
of glassy arrest, as here glassy and fluid dynamics start to de- 
viate. The amplitude h q thus provides important information 
on the physical mechanism causing glassy arrest. It generally 
exhibits a minimum as function of q at the maximum of f q , 
which describes the stucture that gets frozen in at the glass 
transition. 

From simulation data, h q can be obtained via von- 
Schweidler fits to correlators. Table U reveals that in the 
present case these fits require somewhat smaller critical ex- 
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FIG. 5: Enlarged view of the nonergodicity parameters of Fig. [4] 
and structure factors in the low-g-region: Minor differences in f c q 
show up in the region qD < 10 only, where also the input S q differs 
considerably. Simulation data ((D) red symbols) agree qualitatively 
with the model calculation (B) including the barrier. 



portents b compared to all MCT calculations (A-C). This is 
consistent with the observation, see below in Sect. E, that the 
relaxation appears to be somewhat more stretched in the sim- 
ulation than predicted by MCT. 



Von Schweidler's law, %(t) = f c q - h q (t /x) b from Eq. ®, 
gives a much stronger initial relaxation of the frozen-in struc- 
ture for qD < 30 in the simulations than in the MCT calcula- 
tions. In the simulation the small-g modes decay with a larger 
amplitude during the p-process. The overestimate of the sta- 
bility of the incipient glassy structure on length scales larger 
than corresponding to the average particle position indicates 
that MCT misses some of the larger-distance relaxation mech- 
anisms. Nevertheless, the possibility to match a common von 
Schweidler series Eq. ([8]) to the correlators at large and small 
wavevectors 12911 supports our conclusion from Sect. B that 
the structural relaxation for all q is enslaved to local bond for- 
mation. 

The underestimate of h q at low q suggests that MCT un- 
derestimates the contribution of the repulsion driven mech- 
anism of vitrification in the present system. Apparently, in 
the simulatated system attraction driven and repulsion driven 
glass transition compete, and both transition lines are close. 
Our quantitative MCT calculations (erroneously) position the 
system too far from the repulsive transitions. Within MCT, a 
higher order singularity appears in the vicinity of the merging 
of the two glass transition lines, signalled by A, = 1 . Indeed, a 
larger A. exponent is observed in the simulations compared to 
the MCT calculations, which glassify due to bond formation 
solely. 



D) simulation 
A) MCT 
B MCT 
C) MCT 




D. a-relaxation times 



FIG. 6: Critical amplitudes in simulation (D) (red diamonds) and 
MCT (A,B,C) (black dashed-dotted, red full, black dashed) close to 
dynamic arrest. Because of the different oc-times in simulation and 
theory the MCT-results for (A,B) were scaled on the simulation re- 
sults. 

The critical amplitudes in Fig. [6] which are associated with 
the calculations where the attractions dominate (A,B), all ex- 
hibit a very broad peak in q. This shows that very local mo- 
tion takes part in the p-relaxation of an attraction driven glass 
transition. The presence (in (B)) or absence (in (A)) of the pre- 
peak does not influence h q beyond tiny changes for qD < 10. 
MCT thus correctly identifies local bond-formation as the rate 
limiting step during the p-process. MCT underestimates h q in 
the q-range below qD = 30 and thus overestimates the stability 
of the glassy structure on intermediate and long length scales. 



100 
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(B) MCT 
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FIG. 7: a-relaxation times z q in MCT for the systems (A),(B),(C) 
(black dashed-dotted, red full, black dashed line) and in the polydis- 
perse simulation ((D) red diamonds). The curves are normalized to 
the value of z q at qD = 22.7 '. In the inset, the results for models (A) 
and (B) and simulation (D) are shown in a log-log plot. The straight 
line gives the asymptotic behavior %„ = q~^l b valid for large q with 
the corresponding von Schweidler b = 0.37 obtained previously. 



The criterion for quantitatively defining the a-relaxation 
times is somewhat arbitrary. We chose the definition 
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for the a-relaxation times, where in the theory X = 0.05 and 
in the simulation X = 0.50. The latter choice was required 
because of the limited simulation time, but incurs larger cor- 
rections to the values of x q arising from faster relaxation pro- 
cesses. The different definitions are reconciled in the follow- 
ing comparison by normalization of the times, i q /x q o=22.i, 
which brings out the ^-dependence. The a-relaxation times 
x q generally vary in phase with the nonergodicity parameter 
f q and the structure factor S q , a phenomenon often referred 
to as de Gennes narrowing. Repulsion driven glass transitions 
display the largest x q at the principal peak in S q , which in- 
dicates that the cage formed by the particle's next neighbors 
induces the dynamical arrest. Here on the contrary, the slow- 
est relaxation takes place either (i) at the prepeak, when the 
barrier causes void formation, or (ii) for q — ► in model (A) 
without barrier on approaching the phase separation region. 

Figure |7]gives a-relaxation times calculated from the above 
definition and normalized to their value at qD = 22.7 . In MCT 
for models with barrier (B,C) one finds that the slowest modes 
are connected with the prepeak. In model (A) i q decreases by 
more than a factor 2 at this q, when the prepeak in S q is elim- 
inated. The simulations, however, do not allow to check this 
difference, but only show that the slowest modes are those 
with qD < 2. Nevertheless, the different dynamics at small 
q has no further impact on the dynamics at larger q; in the 
inset of Fig. [7] the z q for all models agree at large q, where 
the power-law x q ~ q~ l l b is also tested. It holds nicely in the 
simulation data, with the von Schweidler exponent obtained 
from the fitting of the correlation functions: b = 0.37. MCT 
explains z q for larger wave vectors quite well, though small 
systematic deviations emerge because of the difference in the 
von-Schweidler exponents between simulation and calcula- 
tions, models (A), (B) and (C). Let us note in passing that the 
asymptotic behavior x q ~ q~ l l b holds earlier in the polydis- 
perse MCT calculation (C) (not shown), then in the monodis- 
perse ones (A) and (B), where deviations are still noticeable 
in the inset of Fig. [7] This appears to support the probabilistic 
interpretation of the Kohlrausch law within MCT JMS- 

We conclude that the void structure is completely enslaved 
by the bond formation on local length scales. Even though the 
dissolution of the void structure is the slowest process, there's 
no evidence for a significant influence of the voids on the lo- 
cal dynamics. Local bond formation proceeds identically in 
systems (B,C,D) with barrier and void pre-peak in S q , and in 
system (A) without barrier and void-correlation peak in S q . 



E. Correlation functions 

For an inclusive check of the results beyond asymptotic ex- 
pansions we compare MCT-correlators for a finite distance 
from the critical point obtained as full solutions to the equa- 
tions (fJHU). The fact that in simulations only finite distances 
from the transition point are accessible necessitates this com- 
parison. We used the monodisperse model (B), because it pro- 
vides reasonable nonergodicity parameters and contains the 
barrier like the simulation (D). The wave vectors for the 
MCT-correlators have been chosen to be as close as possible 




FIG. 8: Density correlation functions in simulation ((D) red 
solid lines) and MCT calculations in model (B) (black dashed 
lines): The horizontal black bars indicate the critical plateau val- 
ues of f q in MCT. The wavevectors from bottom to top are qD = 
57.1,45.9,33.9,22.7, 12.5. The short time diffusion coefficient was 
set to D s /D 2 = 0.133. The theoretical curves correspond to a poly- 
mer concentration of § p = 0.2461, which means a separation param- 
eter of about Zmct = —0.001. The simulation results were obtained 
for a polymer concentration fy p = 0.42, which corresponds to a sep- 
aration parameter of £ JI/nli = -0.015. 



to the simulation values (AqD ±0.1). The separation parame- 
ter for the MCT calculation was adjusted so that the structural 
relaxation can be compared most succinctly. 

The correlators from simulation in Fig. [8] show a two-step 
relaxation process with the final a-relaxation from the plateau 
of height f q . The p-process describes the dynamics close to 
f q including some part of the approach to f q from above. The 
simulation data exhibit damped vibrational motion on short 
time scales. This is neglected in Eq. ©, which therefore can 
describe the dynamics only at later times. It is only this struc- 
tural relaxation that MCT addresses and thus the modelling at 
short times is done as simply as possible. We do not attempt 
to (a) include vibrational motion, (b) capture the separation of 
short time and long time dynamics quantitatively, but (c) only 
consider the shape of the structural relaxation in the following 
comparison. 

The results highlight the strong q-dependence of the struc- 
tural relaxation. In addition a strong stretching, i.e. non- 
exponential a-relaxation is also observed. According to dif- 
ferent von-Schweidler exponents of b = 0.37 in the simula- 
tion and b = 0.54 in MCT the simulation results appear some- 
what more stretched. Nevertheless, the local dynamics, where 
the bond formation can be directly seen, is well described by 
MCT. Amplitude and shape of the a-process are rather well 
captured, as would become even clearer if simulations closer 
to the transition could be performed. But for the dynamics on 
larger scales only qualitative statements can be made . 
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VI. CONCLUSIONS 

Glass transitions within MCT are bifurcation points in 
the equations of motion of the structural relaxation. While 
the equilibrium structure of the considered glass forming 
fluid changes smoothly, the dynamics slows down signifi- 
cantly and a metastable glassy structure comes into existence. 
The bifurcation transitions of MCT contain universal signa- 
tures, like von Schweider's law that is the origin of the non- 
exponentiality of the (final or a-) structural relaxation. Non- 
universal amplitudes, like the critical glass form factor f£ or 
the critical amplitude h q entering von Schweidler's law, con- 
tain the information in MCT about the physical mechanisms 
causing arrest. 

We considered the wavector dependence of fj:, h q , and of 
the a relaxation times z q in order to discover the origin of 
glassy arrest in colloidal dispersions of particles with short 
ranged attractions at intermediate packing fractions. The 
study was motivated to gain insight into the connection be- 
tween attraction driven glass transitions at higher densities and 
colloidal gelation at lower ones. We looked at simulations of 
a model system where particles interact additionally with a 
weak long-range repulsive barrier. 

Structure factors S q directly taken from simulations were 
used as the only input to MCT, in order to study the impor- 
tance of the mesoscale gel-like structure indicated by a pre- 
peak in S q . It arises in the simulations from the long-ranged 
repulsive barrier that suppresses gas-liquid demixing. This 
mesoscale peak is the slowest mode in the system, as the void- 
structure takes the longest time to dissolve. Yet, we find that 
the existence of the pre-peak in S q does not affect MCT cal- 



culations for the attraction driven glass transition. The latter 
is caused by local bond formation apparent in all quantities at 
large wavevectors. We conclude that the mesoscopic mode is 
enslaved to the formation of physical bonds, and that this lo- 
cal process is not affected by the larger-scale heterogeneities 
of the system. MCT quantiatively captures local bond for- 
mation but somewhat overestimates the stability of the glass 
on larger length scales. Still, the mechanism of arrest in the 
dispersion at intermediate density is the attraction driven one 
discovered in MCT at higher densities. 

The role of polydispersity in MCT calculations was also 
considered. Averaged structure factors of polydisperse sys- 
tems miss the large g-tail indicative of short ranged attrac- 
tions. This arises from negative interference of the various 
contributions from particle-pairs with different contact dis- 
tances. While this prevents the use of averaged S q from poly- 
disperse systems to capture an attraction driven glass tran- 
sition, it does not imply that the vitrification mechanism in 
polydisperse systems is different. Rather the actual polydis- 
perse system exhibits bond formation, and can be described 
quantitatively within MCT using the appropriate S q from, for 
example, the corresponding monodisperse system. 
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